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Abstract. - We calculate exactly the Casimir force or dispersive force, in the non-retarded 
limit, between a spherical nanoparticle and a substrate beyond the London's or dipolar ap- 
proximation. We find that the force is a non-monotonic function of the distance between the 
sphere and the substrate, such that, it is enhanced by several orders of magnitude as the sphere 
approaches the substrate. Our results do not agree with previous predictions like the Proximity 
theorem approach. 



The technical and experimental advances associated to the investigation of micro and 
nano devices has boosted the research on forces acting at such scale, like Casimir and van der 
Waals (vW) dispersive forces. The origin of dispersive forces between atoms and macroscopic 
bodies may be attributed to electromagnetic interactions between their charge distributions 
induced by quantum vacuum fluctuations, even when they are electrically neutral. In a first 
approximation the charge distribution of neutral particles may be represented by electric 
dipoles. This dipole approximation was employed by London to calculate the non-retarded 
van der Waals interaction potential V v w(z) between two identical polarizable molecules by 
using perturbation theory in quantum mechanics 1 . He found that VVw(z) ~ —a 2 /z e , where 
a is the polarizability of the particles and z is the magnitude of the distance between them. 
In 1948, Casimir and Polder j^j found a correction to the London-vW force by considering 
the influence of retardation, where at large distances the potential varies like 1/z 7 , instead of 
1/z 6 . Furthermore, they found that the vW interaction could be attributed to the change of 
the zero-point energy [3], 

Vvw(z) = £ (z) - £ (z ->• oo) = ^ ^2[uj rn (z) - uj m (z -> oo)], (1) 

m 

where tv m (z) are the classical proper electromagnetic modes of the system. 

Later, Casimir studied a simpler problem |3J: the change of the zero-point energy of 
two parallel conducting plates separated by a distance z, finding an interaction energy per 
unit area, Vc{z) = — (^hc)/ (7 20) (1/z 3 ), where c is the speed of light. In 1956, Lifshitz 
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■4 extended the theory of Casimir to dielectric materials, and later on, it was shown that 
the Lifshitz's formula can be obtained from the zero-point energy of the interacting surface 
plasmons of the plates (SJ|5]. One of the first experimental evidences of the Casimir force 
was observed by Derjaguin and Abrikosova |7|. However, they measured the force between a 
glass plate and spherical lenses to avoid experimental difficulties inherent in keeping the plates 
parallel. To understand their measurements, Derjaguin and collaborators developed the so 
called Proximity Theorem to estimate the Casimir force per unit area between two spheres of 
radii Ri and i?2, in terms of the energy per unit area between parallel plates, Vc(z). Assuming 
that the Casimir force on a small area of one sphere is due to locally "flat" surfaces on the 
other sphere, it is found that 

This theorem holds when z <C Ri,R,2, and it is not clear up to what limit this approach is 
valid. In the limit when i?i = R, and i?2 — > oo, the problem reduces to the case of a sphere 
of radius R and a flat plate, that yields to Fp^(z) = 2nRV(z). If the sphere and plate are 
perfect conductors, 

PT[ ' 360 «3' 

while the lack of retardation implies a behavior of ~ 1/z 2 . 

Casimir and Polder [2] also studied the interaction of a neutral atom with a perfectly 
conducting plate, such that, the force behaves asF~ —a/z 5 at large distances and as —a/z 4 
at small distances (without retardation effects). The same behavior can be also obtained 
calculating the change of the zero-point energy between a induced dipole moment on a sphere 
of polarizability a with its own image dipole in a plate. Contrary to what one may expect, 
this result differs from the one obtained using the Proximity Theorem by a factor of 1/z 2 . 
Therefore, an exact calculation of the dispersive force between a sphere and a planar surface 
becomes essential. In this direction, Langbein [H] used perturbation theory to calculate the 
dispersion forces between two dielectric spheres by considering that they are composed of 
fluctuating electric dipoles. He considered that the dipoles of sphere A interact among them, 
giving rise to a screened field that influenced the dipoles of sphere B, which also induced a 
reaction field on A, and so on. He expanded the field fluctuations in terms of spherical har- 
monics obtaining an infinite Taylor series in terms of R/z. Langbein only evaluated explicitly 
some upper and lower bounds for the interaction energy by considering: (i) few terms of 
the Taylor series, (ii) approximate expressions for the polarizabilties of the spheres, and (iii) 
constant dielectric functions for frequency ranges lower than a cutoff frequency. Under those 
approximations, the upper energy-bound is consistent with the Proximity Theorem. 

More recently, Johansonn and Apell [5] calculated the vW force between a sphere and a 
semi-infinite slab to estimate the errors involved in the Proximity Theorem. They calculated 
the electromagnetic stress tensor in terms of the electric field correlation associated to the 
response to an electric dipole by determining the Green's function of the Poisson's equation 
in bispherical coordinates. They concluded that for small separations the behavior of the 
attractive force is consistent with the Proximity Theorem. However, one drawback in this 
formalism, relevant for the purposes of the present paper, is that in bispherical coordinates 
the section surfaces become planar (a point) for small (large) values of z, so that arbitrary 
values of ratio of z/R cannot be considered. Ford also investigated the Casimir force of a 
dielectric sphere and a wall, within the dipole approximation He found that the relative 
force could oscillate from attractive to negative depending on the distance of the sphere to the 
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wall. However, this oscillatory behavior has not been observed experimentally |Hlll2llT5lll4| 
up to date. 

In this paper, we study the dispersive force between a sphere and a semi-infinite substrate 
in the non-retarded limit. We calculate the interaction energy as the difference of the zero- 
point energy when the sphere of radius R is at a distance z from the substrate, and when 
the sphere is at infinite. The proper electromagnetic modes of the sphere-plate system are 
calculated taking into account all the multipolar excitations, instead of dipolar ones. We 
find these modes using a spectral representation formalism |15| . The spectral representation 
formalism has the advantage that separates the contribution of the dielectric properties of the 
sphere from the contribution of its geometrical properties, such that, it is possible to perform 
a systematic study of the system. 

Our model consist of a nanometric-size sphere of radius R located at a minimum distance 
z from a substrate, with local dielectric functions e sp h(w), and e su b(w), respectively. In the 
non-retarded limit, R and z must be smaller than a characteristic length, I, of the system (for 
a metallic sphere I = c/u) p , where co p is the plasma frequency of the sphere). We assume, like 
in vW, that the fluctuations of vacuum induce a charge distribution on the sphere which also 
induces a charge distribution in the substrate, such that, the Im-th multipolar moment on the 
sphere is given by 

(2/ + 1) j- a( , -rrsubl / O \ 

Qlm = -aim J- \Vi m + V lm ), (2) 

where VJ™ C is the exciting field associated to the quantum vacuum fluctuations at the zero- 
point energy, V^ b is the induced field due to the presence of the substrate, and ai m are the 
polarizabilities of the sphere |16| . Here, the ^m-th multipole moment induced in the sphere is 
defined as 

Qim = [ r' l p(r')Y lm (6',cp')d 3 r', 



where p(r') is the charge density in the sphere, Yi m are the spherical harmonics, and the 
integral is performed over the sphere's volume v s . The ^m-th component of the field can be 
written using a multipolar expansion [T5|, that yields to 

Qlm = aim '^ ^ + E (-ir +l '4m'Ql'm'\, (3) 



4tt \ lm ^ 

V ,m' 

where Qv m ' is the I'm'-th induced multipolar moment in the substrate which is located at 
r = (2(z + R), 9 — ir, (p) from the center of the sphere, and A\™ is the matrix that couples the 
interaction between the multipolar distribution on the sphere and substrate . The induced 
Z'm'-th multipolar moment in the substrate or the "image-multipole" moment is related with 
the Qim on the sphere by 

QVm' ={-l) V+m ' fc{u)Ql-m>, (4) 

where / c (w) is a contrast function that together with (—1)' +m allows to satisfy the boundary 
conditions of the electric field on the plate. The contrast function is given by 

1 + 6 su bM 

and depends only on the dielectric properties of the plate. Substituting eq. Q in eq. (|3J), one 
finds 



E 

I'm' 



4tt Si V 8 mm > , (uj)A i'm' 
[2V + 1) ai' m ' 



Ql'm' = V^ C . (6) 
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The proper electromagnetic modes of the sphere-substrate system satisfy Eq. 10, and are 
independent of the exciting field. Then, the normal frequencies can be obtained when the 
determinant of the matrix in the left-side of Eq. © is equal to zero. One can solve Eq. © by 
a self-consistent procedure, or inverting L(2L + 1) complex matrices of L x L dimension, where 
L is the largest order of the multipolar expansion. To avoid these cumbersome procedures, we 
use instead a spectral representation method to find the proper electromagnetic modes. The 
discussion of the above method can be found in Refll5l 

Suppose that the sphere is homogeneous, such that, its polarizabilities are independent of 
the index to, and are given by 

l[e sph (ui) - 1] 2l+1 

Ctl(UJ) — tt 7 — ; ; K ■ (7) 

' v ' l[e sph {u) + 1] + 1 v ; 

Let us rewrite such polarizabilities like 

= , n \ , R 2l+1 , (8) 

where u(u) = [1 — e sp h(w)] _1 and n l = 1/(21 + 1). The poles of Eq. JSJ, u(u l ) — n l , yield the 
frequencies of the proper electromagnetic modes of the isolated sphere. Therfore, we define 
u(u>) as the spectral variable. Note that in Eq. JSJ), we have separated the material from the 
geometrical properties of the sphere. 

Now, using u(u>), we rewrite Eq. (jjjj as 

V + x w = bp, (9) 
where we have simplified the notation by writing 



Qlm , . (^ 2 ' +1 ) 1/2 T/ 

\i — (/, to), x„= fi , and &„ = — V lm 



vac 



and 



< = 4 v + / e (fl '^ + V ) 1/2 <- (io) 

It was shown in Ref^lthat, is a symmetric matrix that depends on the distance between 
the centers of the"image-multipole" and sphere as, l/[2(z + R)] l+l , such that, H is dimen- 
sionless and symmetric, and depends only on the geometry of the system. Consider the case 
when the contrast factor function f c (u>) is real (we can perform a rotation to the imaginary 
axis (u> — > ioj) to get a real dielectric function), then H is also real. We can always find a 
unitary matrix U that diagonalizes H, such that, 

^(UH)- 1 !!^ =4am v 5 vv ,, 
fin' 

being n v its eigenvalues. The solution of Eq. © is given by x^ = — ^ , G£J b^i, where G£J is 
a Green's operator, whose element can be written in terms of Ujf and n„, as 
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Now, the poles of Gy (it) yield the proper frequencies of the system. Finally, the interaction 
energy is calculated according to Eq. 

The results presented here are calculated as follows. First, we construct and diagonalize 
the matrix H for a given z/R, choosing a maximum value L of multipolar excitations that 
ensures convergence in the interaction energy. Then, we consider an explicit expression for 
the dielectric function of the sphere, and calculate the proper electromagnetic modes trough 



the relation u(uj m ) 



Once we have ui l m , we calculate the energy according with Eq. 



The largest order of the multipolar expansion considered in this work was L = 2000 for 
z/R = 0.001. As a case study, we use the Drude model for the dielectric function of the 
sphere, therefore, u(uj) — [uj(uj + i/r)]/^, where lj p is the plasma frequency and r is the 
relaxation time. We present results for aluminum (Al) spheres with huj p = 15.80 eV, and 
(rw p ) _1 = 0.04. The substrate is sapphire whose dielectric function is real and constant in a 
wide range of the electromagnetic spectrum \17\ . 
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Fig. 1 - Energy as a function of the ratio z/R, when all multipolar effects are taken into account 
(solid line), as well as results when only quadrupolar (dashed line), and dipolar (doted line) effects 
are considered. 

In figm we observe that the energy shows a power law of (z/R)~ 3 when only dipolar 
effects are considered while the curve corresponding to quadrupolar excitations shows three 
different regions: (i) when z/R > 7 the energy shows a power law of (z/R)~ 3 , that means 
that only dipole-dipole interactions are important, (ii) when 2 < z/R < 7 the energy shows a 
power law of (z/R)~ , indicating that dipole-quadrupole interactions are dominant, and (iii) 
when z/R < 2 the energy shows a power law of (z/R)~ 5 , such that quadrupole-quadrupole 
interactions become important. At small distances, the energy of the system can increase up 
to three orders of magnitude when all the contributions from high-multipolar moments are 
considered, as compared with the energy when only dipolar or quadrupolar moments are taken 
into account. Then, the exact energy curve shows a power law of (z/R)~ 3 when z/R > 7, 
when 2 < z/R < 7 the energy shows the same behavior as the curve where quadrupolar effects 
are included, and at smaller distances (z/R < 1) the energy increases sharply as z decreases, 
and it is not possible to assign a power law. This means that the interaction energy is not a 
monotonic function of z like in the Casimir and Polder, and Proximity Theorem approaches. 

In figEl we show the Casimir force for an Al sphere over sapphire calculated as F = 
—d£/dz. We show the force when all multipolar interactions are taken into account, as well 
as up to dipolar, and up to quadrupolar interactions are considered. In agreement with 



6 



EUROPHYSICS LETTERS 




25 50 75 100 

z (nm) 



Fig. 2 - Casimir force as a function of the minimum distance z for a sphere with R — 50 nm. 



the results for the energy, we observed for the force that multipolar effects become evident 
when the minimum distance between the sphere and the substrate is smaller than R. As the 
sphere approaches the substrate, the attractive force suddenly increases up to four orders of 
magnitude as compare with the dipolar interactions. At large distances (z > 2R) the force can 
be obtained exactly if only up to quadrupolar interactions are considered. We also obtained 
that for z > 7R the interaction between the sphere and the substrate can be modeled using 
only the dipolar approximation [15]. 

Shih and Parsegian JHj measured the deflection of atomic beams by gold surfaces to study 
the Casimir-Polder interaction. However, they obtained results that are inconsistent with 
Casimir and Polder, and also with Proximity Theorem, observing larger deflections. Also, 
recent experiments using the sphere-plate configuration have been performed to measure the 
Casimir force ^J^J^HEl an( i their interpretation has relied in the Proximity Theorem. 
However, to compare experimental data to theory, it has been needed to make two significant 
modifications to the Casimir force within the proximity theorem approximation. The first one 
was to employ the Lifshitz formula 0] to calculate the energy density of parallel plants, Vc{z)- 
Then, they employ the Proximity Theorem to calculate the force, Fpt(z). The second modi- 
fication uses the fact that they measured a larger attractive force than the one predicted by 
the Proximity Theorem. They attributed the deviations to the roughness of the surface which 
tend to increase the attraction force like, F(z) = Fpt(z) [l + 6(A r /z) 2 + 15(A r /z) 4 + ■••], 
where the constant A r depends on the model of the surface roughness ^JEII]- ^ n summary, 
the measurements of Shih and Parsegian JS], as well as other experiments ^3], indicate that 
dispersive forces between a polarizable atom or spherical particle and a planar substrate in- 
volves more complicated interactions than the simple dipole model of Casimir and Polder or 
the Proximity Theorem approximation. 

In conclusion, we have shown that the dispersive force between a sphere and a substrate is 
a non-monotonic function of the distance between bodies, showing that the force is enhanced 
by several orders of magnitude as a nanoparticle approaches the plate. On the other hand, 
at large separations the dipolar term dominates the interaction energy, like in the Casimir 
and Polder model. The increment of the force at small separations could explain the physical 
origin of the large deviations observed in the deflection of atomic beams by metallic surfaces, 
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as well as some instabilities detected in micro and nano devices. However, specific experiments 
have to be perform to prove the latter. Our results are in contradiction with the Proximity 
Theorem. As we pointed out before, the Proximity Theorem does not take into account 
the complete geometry of the system. This is more clear if one calculates the zero-point 
energy of the surface modes within the dipolar approximation. In the non-retarded limit, the 
longitudinal surface mode of a dielectric plate is w p i a tc = ^p/V%, while for a sphere the surface 
mode is w sp hcrc = w p / \/3. It is evident that it is not possible to obtain the sphere's surface 
mode from a sum of plate's surface modes, or vice-versa. Then, the proximity theorem is not 
able to include all the effects inherent to the sphere-plate geometry. 

* * * 
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